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=] An abbreviated version of this protocol was published in eLIFE in Jul 2020 
An electrophysiological marker of arousal level in humans 


DOI: 10.7554/eLife.55092 


Detailed protocol 


1. For how to calculate and plot the multitapered power spectrum, refer to the original paper on multitaper spectral analysis: 
Prerau et al. 2016. Sleep trough the lens of multitaper analysis (httos://doi.org/10.1152/physio!.00062.2015) 

They also provide online tutorials and MATLAB code here: 

https://prerau.bwh.harvard.edu/multitaper/ 


2. To calculate the power also see the fieldtrip implementation here: 
https://www.fieldtriptoolbox.org/faq/how_can_i_do_time-frequency analysis_on continuous data/ 
You will find the code in the tutorial. Use "dpss" instead of "hanning". 


3. Code for fits to the PSD is available here 

- poly fit - see code here: https://github.com/voytekresearch/E|Slope 

- robustfit - see code here: https://github.com/jkosciessa/eBOSC 

- model fit (FOOOF) - see code here: https://github.com/fooof-tools/fooof 
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